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ABSTRACT 

We present results from a detailed spectrophotometric analysis of the blue compact dwarf 
galaxy (BCD) Mrk 370, based on deep UBVRI broad-band and Ha narrow-band observations, 
and long-slit and two-dimensional spectroscopy of its brightest knots. 

The spectroscopic data are used to derive the internal extinction, and to compute metallic- 
ities, electronic density and temperature in the knots. By subtracting the contribution of the 
underlying older stellar population, modeled by an exponential function, removing the contri- 
bution from emission lines, and correcting for extinction, we can measure the true colors of the 
young star-forming knots. We show that the colors obtained this way differ significantly from 
those derived without the above corrections, and lead to different estimates of the ages and star- 
forming history of the knots. Using predictions of evolutionary synthesis models, we estimate 
the ages of both the starburst regions and the underlying stellar component. We found that we 
can reproduce the colors of all the knots with an instantaneous burst of star formation and the 
Salpeter initial mass function (imf) with an upper mass limit of 100 M . The resulting ages 
range between 3 and 6 Myrs. The colors of the low surface brightness component are consistent 
with ages larger than 5 Gyr. The kinematic results suggest ordered motion around the major 
axis of the galaxy. 

Subject headings: galaxies: dwarf galaxies - galaxies: starburst - galaxies: compact - galaxies: kinematics 
and dynamics - galaxies: evolution 



1. Introduction 

Blue Compact Dwarf galaxies (BCD) are low- 
luminosity (Mb > — 18 mag), compact (starburst 
diameter < 1 kpc) objects, which have optical 
spectra similar to those presented by H II regions 
in spiral galaxies (Sargent & Searle 1970; Thuan 
& Martin 1981). They are metal deficient galaxies 
— the metallicities of their ionized gas ranging be- 
tween Z Q /50 and Z©/2 — that form stars at high 
rates, able to exhaust their available gas content 
in a time much smaller than the age of the Uni- 
verse. This fact implies that either these galaxies 
are young systems, experiencing their first burst 
of star formation, or that star formation occurs in 
short bursts separated by long quiescent periods 
(Sargent & Searle 1970). Nowadays, the latter ex- 
planation is the more widely accepted, at least for 
the majority of BCDs (Thuan 1991). 

In the recent years BCDs have attracted a 
great deal of interest, and have become key ob- 
jects for understanding fundamental astrophysical 
problems. Much effort has been devoted to the 
study of these galaxies, driven by the following 
fundamental questions: 1) what is the mechanism 
that triggers the starburst activity and is there 
a unique one in all BCD subtypes? 2) what are 
the properties of the underlying stellar population 
(age, metallicity, dynamical status)? Answering 
the above questions would allow to form a gen- 
eral view of dwarf galaxy formation and evolution, 
and in particular to test those evolutionary sce- 
narios linking the different classes of dwarf galax- 
ies, namely, dwarf irregulars (dls), dwarf ellipticals 
(dEs) and BCDs (Thuan 1985; Davies & Phillipps 
1988; Papaderos et al. 1996b; Marlowe et al. 1997; 
Marlowe, Meurer, & Heckman 1999), or the dif- 
ferent morphological classes of BCDs (Noeske et 
al. 2000). However, most of the work carried 
out so far has focused on a statistical analysis 
of BCD samples by means of surface photometry 
(Papaderos et al. 1996a; Telles 1995; Doublier et 
al. 1997; Doublier, Caulet, & Comte 1999; Ostlin 
1998), and has rather demonstrated the complex- 
ity of the topic than given firm answers. 

Prompted by this lack of conclusive results, we 
started an extensive study of a large sample of 
BCDs. The first step was the analysis of broad- 
band observations in B, V , R and / for a sam- 
ple of 28 BCDs. In Cairos et al. (2001a, hereafter 



Paper I) we presented deep contour maps in the 
B band, and surface brightness and color profiles, 
which were used to examine and discuss their mor- 
phology and structure. In Cairos et al. (2001b, 
hereafter Paper II) we enlarged the dataset by in- 
cluding U and Ha observations; we provided inte- 
grated photometry of the galaxies and produced 
an atlas of detailed color and Ha maps. The qual- 
ity of the data allowed us to identify the different 
star-forming regions, and discriminate them from 
the regions where the star formation has already 
ceased. 

One of our main conclusions was that opti- 
cal broad-band photometry alone does not allow 
to disentangle the distinct stellar components in 
BCDs and derive their evolutionary status. The 
measured broad-band colors are strongly affected 
by interstellar reddening and gaseous emission; be- 
sides, we need to separate the light coming from 
young stars from the emission of those born in pre- 
vious star-formation episodes — a methodological 
aspect often overlooked. Spectroscopic informa- 
tion is required in order to correct for interstellar 
reddening and emission line contamination, and 
JZ-band and Ha measurements are fundamental, 
since they are much more sensitive to the star- 
formation properties. Nevertheless, such kind of 
analysis is scarce in the literature, mainly because 
it requires a lot of observing time. In particular, 
not enough effort has been devoted so far to the 
characterization of the stellar population hosting 
the actual burst, indeed a challenging task. 

For these reasons, we decided to perform a more 
detailed study of a number of individual galax- 
ies, based on a new strategy that combines broad- 
and narrow-band photometry with long-slit and 2- 
D spectroscopy. The objects were selected among 
those presented in Papers I and II, so as to form a 
sample including the different subtypes of objects 
classified as BCDs — a not so well defined class 
which spans a wide range of luminosities, metal- 
licities and morphologies (Kunth & Ostlin 2000; 
Cairos 2000; Cairos et al. 2001b). 

This paper is the first in a series which present 
the results obtained from the use of this technique 
to analyze the selected galaxies. Here we describe 
the technique and its application to the galaxy 
Mrk 370, a luminous BCD, with M B = -17.20 
(Paper II), located at a distance of 12.9 Mpc 
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(Thuan & Martin 1981) 1 . 

We have chosen this galaxy in order to study 
a prototype of the group of BCDs populating the 
higher luminosity and metallicity range within the 
most common BCD morphological subclass, the iE 
type (Loose & Thuan 1986, hereafter LT86). This 
type is characterized by an underlying, low bright- 
ness elliptical component which hosts numerous 
star-forming regions. 

In Mrk 370, the brighter knots of star- formation 
are located in the inner part of the galaxy, from 
where smaller star-forming knots emerge and form 
a structure resembling spiral arms (see Papers I 
and II) . The choice of Mrk 370 rests also on the 
fact that the comparison of the properties of its 
many individual starburst knots provides valuable 
information on the mechanisms that trigger the 
star formation, allowing to study how it prop- 
agates within the galaxy, and is a possible test 
of the stochastic self-propagating star formation 
mechanism (Gerola, Seiden, & Schulman 1980), 
which has been suggested to explain how the star- 
formation process in dwarfs is maintained (Coziol 
1996). 

In addition, Mrk 370 may be seen as a pro- 
totype of the BCDs classified as multiple nu- 
cleus starburst galaxies. The presence of two cen- 
tral star-forming knots led Mazzarella & Boroson 
(1993) to label this galaxy as a "double nucleus" 
BCD, prompting further suggestions of its merger 
origin. As we will show in this work, the kinematic 
and spectrophotometric results obtained cast seri- 
ous doubts on the merger scenario. 

A complete collection of color and Ha maps 
of this galaxy can be found at the following Web 
page: http://www.iac.es/proyect/GEFE/BCDs. 
The basic data of Mrk 370 are shown in Table 1. 

2. Observations and Data Reduction 
2.1. Long-slit spectra 

Long-slit spectroscopy of Mrk 370 was obtained 
in September 1998, with the Isaac Newton Tele- 



For consistency with Papers I and II, we use here the dis- 
tance from Thuan & Martin (1981), which is based on 
the radial velocity with respect to the Local Group and 
Ho = 75 km sec -1 Mpc^ 1 ; applying the correction for 
LG infall into Virgo would give D = 10.9 Mpc, and total 
luminosities, absolute magnitudes and sizes would change 
accordingly. 



scope (INT) at the Observatorio del Roque de Los 
Muchachos (ORM, La Palma). The data were 
taken using the Intermediate Dispersion Spectro- 
graph (IDS), equipped with a CCD. A slit 3'8 long 
and l'/2 wide was used. The position angle of the 
slit was set to 125°, so as to pass through the two 
brightest central knots. The spectral sampling was 
1.85 A per pixel, the wavelength range 3600 — 7300 
A, and the spectral resolution ~ 5.5 A (FWHM); 
the spatial sampling was 0.39 arcsec per pixel. The 
seeing was about l'/5. The log of the observation 
is shown in Table 2. 

Data reduction was performed using IRAF 
standard tasks: after subtracting the bias and 
flat-fielding the spectra, they were calibrated in 
wavelength. The sky spectrum was derived by 
averaging the signal in two windows 30-40 pixel 
wide outside the region where the object is still 
detectable, and subtracted out. The spectra were 
then corrected for atmospheric extinction and cal- 
ibrated in flux by means of observations of spec- 
trophotometric standards. 

2.2. Two-dimensional spectroscopy 

The 4.2m William Herschel Telescope (WHT, 
at the ORM, La Palma), equipped with the fiber 
system INTEGRAL (Arribas et al. 1998) and the 
WYFFOS spectrograph (Bingham et al. 1994), 
was used to observe the central region of Mrk 370. 

INTEGRAL currently has three fiber bundles 
with different spatial configurations on the focal 
plane. These three bundles can be interchanged 
during the observations depending on the scien- 
tific program or the seeing conditions. On the fo- 
cal plane, the fibers of each bundle are arranged in 
two groups: one forms a rectangle which maps the 
object, the other forms an outer ring which sam- 
ples the sky-background. Mrk 370 observations 
were done with the standard bundle 2 (hereafter 
SB2). It consists of 219 fibers, each Of! 9 in diam- 
eter; the central rectangle is formed by 189 fibers 
covering an area of 16" x 12'.'3, while the other 
30 fibers are evenly spaced on a ring 90" in di- 
ameter. Figure 1 shows the actual distribution of 
those fibers in the focal plane. At the entrance of 
the WYFFOS spectrograph, the fibers are aligned 
so as to form a pseudo-slit. 

We observed Mrk 370 on August 22 1999, with 
a seeing of ~ 1"3. The spectrograph was equipped 
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with a grating of 600 groove/mm, and a Tek6 CCD 
array of 1124x1124 24-micron pixels. The spectral 
range coverage was 4400 — 7150 A, with a spectral 
resolution of about 6 A; it includes the H/3, [O III], 
Ha, [N II], and [S II] emission lines (see Appendix 
A). 

Data reduction was performed within the IRAF 
environment, and includes bias and scattered light 
subtraction, fiber extraction, fiber throughput cor- 
rection, wavelength calibration, rejection of cosmic 
ray events, and correction from differential atmo- 
spheric refraction effects (for details on reduction 
of fiber data see e.g. Garcfa-Lorenzo, Mediavilla, 
& Arribas 1999, and references therein). 

In Appendix A we present the nuclear spectrum 
of Mrk 370 in the full wavelength range observed 
with the fiber system. We also show the individ- 
ual spectra corresponding to each of the 189 ob- 
served positions (fibers) in a selected spectral in- 
terval which includes the most important emission 
lines. 

2.3. Ha imaging 

Narrow-band images centered on the Ha line 
and on the adjacent continuum were taken in 
1997 December at the 2.2-m telescope of the 
German-Spanish Astronomical Observatory on 
Calar Alto (Almcria, Spain). The instrumentation 
consisted of the Calar Alto Faint Object Spectro- 
graph, CAFOS, and a 2048 x 2048 SiTe CCD chip. 
With a pixel size of 0.53 arcsec, it provides a field 
16 arcmin wide; however, due to the physical size 
of the available filters, only the central round area 
of about 11 arcmin is free from vignetting. The 
average seeing was 1.5 arcsec. 

The image reduction was carried out using 
IRAF. Each image was corrected for bias, using 
an average bias frame, and was flattened by divid- 
ing by a mean twilight flatficld image. The frames 
were then registered (for each filter we took a set of 
dithered exposures) and combined to obtain the fi- 
nal frame, with cosmic ray events removed and bad 
pixels cleaned out. The average sky level was esti- 
mated by computing the mean value within several 
boxes surrounding the object, and subtracted out 
as a constant. Flux calibration was done through 
the observation of spectrophotometric stars from 
the lists in Oke (1990). 

The complete log of the observations is shown 



in Table 2. 

2.4. Broadband Imaging 

UBVRI imaging of Mrk 370 was carried out in 
October 1998 at the 1-m Jacobus Kapteyn Tele- 
scope, JKT (ORM, La Palma). We collected CCD 
images at the f/15 Cassegrain focus, using a Tek 
1024 x 1024 pixel chip, yielding a pixel size of 0.33 
arcsec and a field of view of 5.6 x 5.6 arcmin. The 
average seeing was l'/5. 

A log of the observations is given in Table 2. 
The image processing was performed using stan- 
dard procedures available in IRAF. Flux calibra- 
tion was done through the observation of photo- 
metric stars from Landolt (1992) throughout the 
night. 

3. Spectroscopic Analysis 
3.1. Long-slit spectroscopy 

The position of the slit we used is shown in 
Figure 2, overplotted on a contour map of the 
continuum-subtracted Ha image. (The meaning 
of the labels is explained in Section 4). The spatial 
run of the intensity of the brightest emission lines 
(Ha, H/3, [O III]) has been used as a base for defin- 
ing the regions from where the one-dimensional 
spectra of each knot were extracted by integrating 
the signal within them. Two distinct starburst 
regions have been identified, marked A and B in 
Figure 2. 

Figure 3 shows the summed spectra of the two 
knots, and the spectrum integrated over the whole 
galaxy (which we shall refer to as s). The spec- 
trum of region A shows a high continuum level. 
Absorption wings are resolved in H/3 and H7, 
H<5 being visible in absorption only. The Balmcr 
discontinuity, together with these absorption fea- 
tures, can be interpreted as the product of a sub- 
stantial population of older stars, with ages con- 
siderably larger than 10 million years. Region B 
shows an almost flat spectrum, characteristic of 
a dominant OB population, with no evident ab- 
sorption features. The integrated spectrum, S, is 
similar to that of region A. 

Fluxes and equivalent widths of the emission 
lines were measured using the gaussian profile fit- 
ting option in the iraf task SPLOT (a direct inte- 
gration of the flux inside each line gave very simi- 
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lar results). It is known that these measurements 
are an underestimate of the real flux of the Balmer 
lines, because of the underlying absorption compo- 
nent. In the present spectra we could not deblcnd 
the two components. To correct for the underlying 
stellar absorption, some authors (McCall, Rybski 
& Shields 1985; Skillman & Kennicutt 1993) adopt 
a constant equivalent width (1.5 — 2 A) for all the 
hydrogen absorption lines. However, the actual 
value of their equivalent width is uncertain and 
depends on the age of the star formation burst 
(Diaz 1988; Cananzi, Augarde, & Lequeux 1993; 
Olofsson 1995). Thus we preferred to leave it as a 
free parameter, which was determined as follows. 

We assume that the equivalent width, W a b s , of 
the absorption line is the same for all the Balmer 
lines, but can vary from knot to knot. We then 
start with an initial guess for Wab s , correct the 
measured fluxes, and determine the extinction co- 
efficient at A = 4861 A, C(H/3), and its uncer- 
tainty, AC (H/3), through a fit to the Balmer decre- 
ment, given by the equation: 



F(X) /(A) 



F{Hp) I {Hp) 

F(X) 



x 10 C(ff/3)x/(A) 



(1) 



where is the line flux corrected for absorp- 

tion and normalized to H/3; jr^m is the theoretical 
value for case B recombination, from Brocklehurst 
(1971), and /(A) is the reddening curve normal- 
ized to H/3 which we took from Whitford (1958). 

We then search, for each knot, the value of Wabs 
that provides the best match — that is, the lowest 
AC (H/3) — between the corrected and the theo- 
retical line ratios. 

Reddening-corrected intensity ratios and equiv- 
alent widths are quoted in Table 3. 

The [O III] A 4363 line could not be detected 
in the spectra, so the oxygen abundance for the 
different regions of Mrk 370 was derived using the 
empirical 12+log O/H vs R23 relation, first de- 
rived by Pagel et al. (1979). We have used the R23 
calibration according to McGaugh (1991). The 
ratio [N II] /[O II] was used to break the degen- 
eracy in the calibration in order to select either 
the upper or the lower branch (McGaugh 1991). 
Since the values of [N II] /[O II] for the three spec- 
tra analyzed are always above the 0.1 threshold 
(logA6584/A3727 = -0.59, -0.81, -0.75 for re- 



gions A, B and s), the upper branch of the calibra- 
tion has been used for the determination of the 
abundance. We have derived the following values: 
12 + \ogO/H = 8.7, 8.6 and 8.7 for the spectra 
of regions A, B and S respectively, for which the 
corresponding values for the logarithm of electron 
density are 2.07, 2.2 and 2.3. 

The average value for the nitrogen to oxygen 
ratio is log N/O = —1.1. The derived average 
oxygen abundance is approximately Zq/2. The 
abundances are equal for regions A and B, within 
the uncertainties. 

3.2. Two-dimensional spectroscopy 

Two-dimensional spectroscopy with optical 
fibers allows to collect simultaneously the spectra 
of many different regions of an extended object, 
combining photometry and spectroscopy in the 
same dataset. 

From the individual measurements (line fluxes 
and ratios, continuum level, radial velocity, etc.) 
obtained at each fiber, we can build 2-D maps 
by using a two-dimensional interpolation tech- 
nique, such as the Renka & Cline algorithm (im- 
plemented in the E01SAF and E01SBF routines 
in the NAG Fortran Library). This requires a pre- 
cise knowledge of the exact location of each fiber, 
which is determined by a metrology machine. In 
this way we built up images of 95 x 95 pixels with 
a scale ~ Of! 2 pixel -1 . While the spatial sam- 
pling (that is, the fiber diameter) of INTEGRAL 
(STD2) is 0.9 arcsec, the centroid of any peak in 
our maps can be measured with an accuracy of 
around 0.2 arcsec (that is 1/5 of the fiber diame- 
ter; see for instance Mediavilla et al. 1998) 

3.2.1. Continuum and Line- Intensity Maps 

Because of non-photometric conditions dur- 
ing INTEGRAL observations, we could not flux- 
calibrate directly our spectra. The flux calibration 
was done by using the long-slit data as a refer- 
ence. We integrated the signal in our 2-D spectra 
within two regions matching regions A and B in the 
long-slit spectra. Comparing the fluxes of differ- 
ent lines, we could compute the conversion factor 
with an accuracy of about 20%. This is the error 
associated to the fluxes derived from INTEGRAL 
data using this calibration. 

Figure 4(a) shows a continuum map (within 
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a spectral region where no emission lines are 
present) of Mrk 370 from INTEGRAL spectra. 
We also present in Figure 4(b) a V band image 
of the galaxy, in order to evidence the excellent 
agreement between the continuum map and the 
broad-band images. 

The INTEGRAL data show that the position of 
the peak in the continuum maps does not depend 
on what spectral range was used to build them. 
Therefore, we can safely define this point as the 
" optical nucleus" of the galaxy, and we set it as the 
origin of the coordinate system in all our spectral 
maps. 

By fitting a single-gaussian to each emission line 
in our spectra, we have obtained emission line flux 
maps. Figure 4 (c) shows the Ha emission map. 
We also fitted an absorption component to Balmcr 
lines in those spatial positions where the absorp- 
tion wings were clear. For comparison, an Ha nar- 
row band filter image of Mrk 370 is displayed in 
Figure 4 (d). The area covered by INTEGRAL in- 
cludes four of the many star-forming knots present 
in the galaxy. The optical nucleus does not coin- 
cide with any of them. 

Figure 4 (e) is the map of the [O III]/H/3 ratio, 
which is commonly used as an indicator of the ex- 
citation. The morphology of the [O III]/H/3 map is 
similar to the emission lines intensity maps. The 
regions with the highest ionization coincide with 
the knot peaks, and the level of excitation is nor- 
mal for star-forming regions. Figure 4 (f) is the 
map of the extinction coefficient C (H/3) computed 
from the observed Ha/H/3 ratio. 

The peaks in the extinction coefficient map are 
close to, but not coincident with, the peaks of the 
Ha knots. The average offset for the four knots in 
our observed area is 0.26 ±0.11 arcsec. While this 
displacement is small — in fact, it is at our de- 
tection limit — , it may suggest a real spatial shift 
between the center of the HII regions and the holes 
in the C(Ha) distribution. A similar case has been 
reported before by Mafz-Apellaniz et al. (1998) for 
the starburst knots of the galaxy NGC 4214. They 
find the dust to be concentrated at the boundaries 
of the ionized regions. The explanation that they 
propose is that stellar winds have contributed to 
diffuse the dust into the inter-cluster medium. 



3.2.2. Kinematical Pattern 

The kinematic analysis is limited to the emis- 
sion features, as absorption lines are too weak to 
be of some utility. 
Individual Radial Velocities 

Radial velocities were obtained by fitting a sin- 
gle gaussian to each emission line present in our 
spectra. Although the spectral resolution of the 
Integral Field Spectroscopy is ~ 275 km/s, the 
centroid of the emission features could be mea- 
sured with an accuracy of ±10 km/s, as indicated 
by the scatter among the radial velocities obtained 
from different lines. The final radial velocities of 
the ionized gas were obtained by averaging the re- 
sults derived from individual emission lines. We 
compared the results for the two different indi- 
vidual exposures of Mrk 370. The mean differ- 
ence among the individual values was 0± 12 km/s. 
Therefore, by adding the error on the wavelength 
calibration (<~ 7 km/s), we estimate that the fi- 
nal uncertainties on the velocity measurements are 
< 20 km/s. 

The Mean Ionized Gas Velocity Field 

Figure 5 shows that the velocity field appears to 
be regular and similar to that of a rotating disk. 
The west region is receding, the east part is ap- 
proaching. In order to determine the location of 
the kinematic center (k hereafter), we computed 
the "differential" of the velocity field, averaging 
the absolute differences between one point and its 
nearest neighbors. This procedure emphasizes re- 
gions in which the velocity is changing rapidly. 
The position derived in this way is placed ~ 2'.'4 
SE (-2"36, Of.' 37) of the optical nucleus. K is lo- 
cated between regions A and B, in an apparently 
non-significant location in the galaxy, and has a 
velocity value of 766±20 km/s, consistent with the 
systemic velocity of the galaxy, 783±16 km/s, re- 
ported by Falco et al. (2000). 

K is nearly coinciding with the geometrical cen- 
ter of the external isophotes, obtained from our 
broad band images. This fact might suggest that 
the star-forming knots in the central region of 
Mrk 370 could share the same velocity field of the 
main body. 

Velocity Field Parameters 

To quantify the relevant parameters of the ion- 
ized gas velocity field, we adopted a simple kinc- 
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matic model (Mihalas & Binney 1981): 

V T (R,0) = V sys +Z(R)cosi+Q{R)cos[6-5(R)]siTii 

(2) 

where V sys is the systemic velocity Q(R) is the 
strength of the velocity field projection onto the 
galactic plane, and 6(R) is the PA of the major 
kinematic axis. R and 9 are the galactocentric 
polar coordinates in the plane of the galaxy. 

We fit this model to the data, keeping fixed the 
position of the center at K, and assuming an initial 
estimate of the position angle of the kinematical 
major axis So ~ 90°, and an inclination of the 
galaxy i <~ 40° (Nordgren et al. 1995). The results 
do not depend on So, and are rather insensitive to 
reasonable changes in i. In Figure 6 we present 
the dependence on galactocentric distance of the 
parameters resulting from the fit. The systemic 
velocity we derived is 761 km/s, with an internal 
error of ±2 km/s. 

4. Photometric Analysis 

4.1. Ha photometry of the starburst knots 

The left panel of Figure 7 shows a grey-scale 
Ha image of Mrk 370. The galaxy reveals a com- 
plex structure consisting of several resolved knots 
located atop an extended low surface brightness 
(LSB) component having regular appearance (see 
the R band map in Figure 7, right panel) and sub- 
stantially redder colors (Paper I and Paper II). 

In the central part of the galaxy are two large 
star-forming regions (labeled 11 and 12 in Fig- 
ure 2), « 7 arcsec (450 pc) apart, which appear to 
be connected by a faint bridge- like structure; these 
two star-forming regions have been previously cat- 
aloged as " double nucleus" by Mazzarella & Boro- 
son (1993). Smaller star-forming knots extend out 
in the north-east and south-west directions, to ap- 
proximately 40 arcsec (2.5 kpc) from the optical 
center (i.e., the center of the outer isophotes). The 
star-formation activity is concentrated along this 
direction, with the exception of region 11). 

We have used the FOCAS package to identify 
and analyze the different knots in the continuum- 
subtracted Ha images. The FOCAS task looks for 
local maxima and minima in the ADU counts of 
each pixel in the image. We have adopted the 
criteria that, for a knot to be detected, its area 
must be larger than 12 pixels (corresponding to 



an equivalent diameter of 2 arcsec or 125 pc), to 
ensure that the diameter of the knot is larger than 
the PSF, and the counts of each of its pixels must 
be higher than 3 times the standard deviation of 
the sky background. Using these criteria, a to- 
tal of 16 knots were detected; they are labeled in 
Figure 2. 2 

Data on the Ha flux, luminosity, number of ion- 
izing photons, ionized hydrogen mass and equiv- 
alent width of each knot were calculated; these 
results are presented in Table 4. The number of 
ionizing photons and the ionized hydrogen mass 
were calculated following Kennicutt (1988). The 
electron temperature and density needed to calcu- 
late the mass of the hydrogen ionized gas and the 
number of ionizing photons were obtained from 
the long-slit spectra. 

The resulting parameters have been corrected 
for interstellar extinction and for the contribution 
of the [N II] lines. For the Knots 11 and 12, we 
used the spectroscopic parameters derived for the 
subregions B and A (see Figure 2). We used the 
bi-dimcnsional data to correct knots inside the in- 
tegral field. For those knots outside the region 
covered by our 2-D spectra, the above corrections 
have been applied by using averaged values de- 
rived from the central region. 

The luminosity of the brightest knots (11 and 
12) are comparable to those of individual giant 
H II regions (typically 10 39 -10 40 erg sec" 1 ; Ken- 
nicutt 1991). The range in luminosities is similar 
to that found for the H II regions cataloged in the 
Magellanic Irregular NGC 4449 (Fuentes-Masip et 
al. 2000). 

The equivalent width VF(Ha) was computed as: 

VF(Ha) = C Ra /C cont x AWjjq, (3) 

where CjJq. is the total number of net counts of 
the knot, C cont is the corresponding total number 
of counts in the rescaled continuum image, and 
AWjjq, is the width of the Ha filter transmission 
curve. 

The Ha equivalent width has been corrected for 
the contribution of the underlying stellar emission 
to the continuum level; it means that C cont are 
the total counts measured in the continuum band 



2 The Integral continuum maps show another fainter knot, 
not detected by Focas as it does not satisfy these criteria. 
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after subtracting the counts due to the emission 
of the older stars. In order to estimate this latter 
contribution, we assumed that the underlying pop- 
ulation can be described by an exponential model 
(this procedure is explained in the next section). 
As the exposures in the continuum narrow-band 
filter were not deep enough to derive the properties 
of the underlying disk, the structural parameters 
derived in R band have been used. The continuum 
flux map through the Ha filter has been derived 
from the R frame after the latter has been rescaled 
so as to match the continuum in the former. We 
followed the same procedure, based on field stars, 
commonly employed to match Ha and its adjacent 
continuum (see Paper II for a description of this 
technique) . 

In Table 4 we show the values of W(Ha) of the 
different star-forming regions in Mrk 370; we have 
listed the values before and after applying the cor- 
rection for the contribution of the older stars, to 
provide a direct estimation of the importance of 
such contribution. We found the difference to be 
significant in most cases. 

4.2. Broad-band Photometry 

4.2.1. The Starburst 

The young starburst regions have been isolated 
in each broad-band filter image. Using FOCAS, we 
built a mask where the emission regions delimited 
in the Ha frame have the constant value one and 
the rest are set to zero. Each broad-band image is 
then multiplied by this mask. Finally we perform 
aperture photometry on these images, to derive 
magnitudes and colors of each knot. 

Table 5 shows the broad-band photometry of 
the sixteen knots detected in Mrk 370; the mea- 
surements have been corrected from Galactic ex- 
tinction using the data published by Burstein & 
Heiles (1984). The extinction coefficient in the 
other bands were related to the B value following 
Rieke & Lcbofsky (1985). 

In order to characterize the stellar population 
of the individual regions, we compare their inte- 
grated properties with the predictions of evolu- 
tionary synthesis models. This task is especially 
difficult in this case, because at least three differ- 
ent contributors are producing the observed light: 
the young stars, the gaseous nebula surrounding 
them, and the older stars underlying the ongoing 



burst. Broad-band colors must be corrected from 
the contribution of these different factors before 
they can be compared with the model predictions. 
In order to derive the "true" colors of the star- 
burst, we proceeded in the following way. 

First, colors were corrected for interstellar red- 
dening, using the value of C(H/3) derived from our 
long slit spectra. For those knots lacking spectro- 
scopic data, we used the C(H/3) averaged over the 
whole long-slit spectrum (which is probably our 
best possible guess). 

Second, we have corrected for the contribution 
of emission lines. In gas rich objects, like BCDs, 
the presence of strong emission lines in the spec- 
tra can severely affect broad-band colors; the exact 
amount of the contribution of the lines to the flux 
observed in a given filter depends on the equivalent 
width of the line and on its precise location under 
the filter transmission profile. We used our spec- 
troscopic information to compute the fraction of 
flux in the observed broad-band filters accounted 
for by each emission line, by using the task cal- 
cphot in the synphot package of STSDAS. 

Last, we subtracted the emission of the under- 
lying population of stars. We assumed that the 
galaxy is composed of the starburst regions delim- 
ited by the FOCAS task, and of an underlying low 
surface brightness component, which is described 
by an exponential function (see next section). 

The integrated magnitude of a knot, in a given 
band, A, is given by: 

m A = -2.5 x log (Cf B + Cf ) + K\ (4) 

where Cf B is the number of counts in the knot 
coming from the starburst; is the contribu- 
tion of the underlying disk within the knot region, 
and K\ is the corresponding calibration constant 
(which incorporates exposure time, airmass, etc.). 

is easily computed from the central surface 
brightness and scale length of the best-fit expo- 
nential model, and subtracted out in the previous 
formula, yielding the corrected magnitude of each 
knot: 

m x (SB) = -2.5 x log (Cf B ) + K x (5) 

This process has been applied to the B, V, R 
and / bands. The images in U were not deep 
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enough to allow a reliable fit to the light profiles; 
however, by using typical values oiU—B for an 
evolved stellar population (Bruzual 1993; Vazdekis 
et al. 1996), we have estimated the contribution of 
the host galaxy in the U band. We found that, in 
this pass-band, the older component accounts for 
only ~ 1% of the knot flux, so we can neglect 
this correction. The resulting corrected colors are 
listed in Table 5. 

We would like to stress that uncorrected and 
corrected knot colors differ significantly. In the 
case of the Mrk 370, while emission lines do not 
change much the broad band colors (shifts in 
B — V arc w 0.05 or less), the internal extinction 
may change B — V by up to 0.25 mag, and the 
effect of the contribution of the LSB stellar com- 
ponent is equally important. The final color shift 
is clearly function of the position in the galaxy; 
this means that it is dangerous to interpret color 
gradients as genuine variations of the knot proper- 
ties (for instance their age) if the above corrections 
have not been applied. 

4-2.2. The Underlying Stellar Population 

Surface brightness profiles (SBP) in BVRI, as 
well as color profiles of Mrk 370, have been pre- 
sented in Paper I. The SBP of the galaxy was used 
to decompose the galaxy in two components: the 
starburst, Ha emitting regions, plus an older LSB 
stellar host. By fitting an exponential function to 
the faint, outer parts of the galaxy light profile 
(where we expect no or only negligible contribu- 
tion from the starburst), we determined the struc- 
tural parameters (scale length a and central sur- 
face brightness /t*o) of the LSB host. These param- 
eters were used to calculate its total magnitudes 
in the different bands. 

The structural parameters in the different 
bands, together with the total luminosities of the 
starburst and of the old stellar component, are 
given in Table 6. 

Colors of the underlying component were de- 
rived from the best-fit exponential model; we 
found {B-V) host = 0.80, (V-R) host = 0.50 and 
(V-I) host = 1.46. 



5. Evolutionary properties of Mrk 370 
5.1. The Starburst 

In order to constrain the properties of the indi- 
vidual H II regions, we compare their observed 
properties with the predictions of evolutionary 
population synthesis models. These models assume 
that stars are born with masses distributed ac- 
cording to the initial mass function (imf), with a 
star formation rate (sfr). The stars are evolved 
according to theoretical evolutionary tracks and 
finally, using empirical or theoretical calibrations, 
the predicted colors or spectra of the composite 
population are obtained. In the last years, peo- 
ple have been doing a lot of work in evolutionary 
synthesis models, see for instance Leitherer et al. 
(1996). Depending on the hypothesis on which 
the different models are based, the predicted val- 
ues for the observables can be slightly different (for 
example, theoretical modeling of the evolution of 
massive star is still undergoing significant changes, 
Maeder & Conti 1994; and revision of these models 
often leads to new generations of synthesis calcu- 
lations) . This means that it is important to exam- 
ine the characteristics of the available models and 
select the one most suitable to each specific case. 

We have compared the derived observables of 
the starburst with the predictions of Leitherer et 
al. (1999, hereafter SB99) evolutionary synthesis 
models, known as Starburst99. These models pro- 
vide a wide set of predictions for the spectropho- 
tometric properties of star-forming galaxies. They 
are an improved version of the ones previously 
published by Leitherer & Hcckman (1995): the 
latest set implements the stellar evolution mod- 
els of the Geneva group, as well as the model at- 
mosphere grid compiled by Lejeune, Cuisinier, & 
Buser (1997). 

For each of the sixteen delimited knots, we 
searched for the model which best matches the ob- 
served parameters: W(Ha), U — B, B — V, V — R 
and V—I. In order to break the age and metal- 
licity degeneration, we set the metallicity equal to 
that derived from the emission line fluxes (which 
is a reasonable approximation to the metallicity of 
a young population) . 

We found that, for all the knots, we can re- 
produce the observed quantities with an instanta- 
neous burst (IB) of star formation and the Salpeter 
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IMF with an upper mass limit of 100 Mq. No 
restriction can be put on the lower mass, given 
its poor sensitivity to the young burst evolution, 
which is dominated by the contribution of massive 
stars (M > 8M ). The derived ages from U-B 
and Ha are quoted in Table 7; both observables 
give usually ages consistent with each other. 

The ages derived using B — V, V — R and V — I 
are slightly larger, but in general compatible with 
the former. We notice that these colors are not 
very sensitive to the evolution of the young stars, 
as they vary by about 0.3 to 0.4 mag when the 
age increases from to 100 Myr (SB99). Also, 
these colors are strongly contaminated by the dif- 
ferent factors described in Section 4.2.1 (redden- 
ing, contribution of emission lines and emission of 
the older stars). 

Although we have attempted to correct for such 
factors, the results might be imprecise due to 
the various hypotheses we made to work them 
out. First, we adopted a constant C(H/3) in those 
zone having no spectral information, while there 
may be spatial variations in the dust distribution, 
and C(H/3) can vary from knot to knot. That 
this is likely the case is shown by the results we 
obtained for the two central knots: C(H/3)a = 
0.31 and C(H/3)b = 0.39, which translates into 
E(B - V) A = 0.22, and E(B - V) B = 0.27, and 
E(V - I) A = 0.34 and E(V - I) B = 0.43. If 
such variations are detected in two knots so close 
to each other, we can expect larger differences for 
the more peripheral knots. The same result is ob- 
tained from 2D spectroscopy: the 2D distribution 
of C(H/3) do not present an homogeneous mor- 
phology (see Figure 4), but larger C(H/3) close to 
the center of the star-forming knots. 

The same applies to the contribution of emis- 
sion lines, though, as mentioned earlier, they play 
a secondary role. 

Furthermore, the uncertainties in the structural 
parameters obtained by the exponential fit of the 
LSB component also translate into uncertainties 
in the resultant colors. 

All in all, the combined errors on the knot col- 
ors due to the various factors here discussed can be 
easily as large as 0.2 — 0.3 mag; this fact, together 
with their low sensitivity to the starburst proper- 
ties, imply that B — V, V — R and V — I alone can 
not constrain the properties of the young stellar 



population. 

On the other hand, U — B and Ha are much 
more sensitive to the evolution of the young stars. 
They drastically change during the first million 
years of the stars life (see Figures 55, 56, 83 and 
84 in SB99), which converts them in reliable indi- 
cators of the properties of star- forming regions. 

We remark that none of the observable param- 
eters will constrain accurately the properties of a 
star-forming region, unless we correct properly for 
the factors described above. The parameters mea- 
sured directly (i.e. without correction) usually do 
not fit the models, and, if they do, they tend to 
give much older ages. This point is illustrated in 
Figures 8 and 9. 

5.2. Age of the older stars 

The V—I color appears to be much too red to 
be compatible with a normal older stellar popula- 
tion. The data in the / band are not as good and 
deep as in the other bands, and the large uncer- 
tainties on the /-band structural parameters may 
be responsible for this discrepancy. Thus, the fi- 
nal ages have been derived by using the B — V and 
V — R colors. These colors have not been corrected 
from interstellar extinction, because the extinction 
coefficients we derived were based on the gas emis- 
sion lines, and they may not apply to the regions 
outside the area occupied by the starburst. We 
have assumed that the metallicity of the stellar 
population is lower than the metallicity derived 
for the gas. 

We have compared the B — V and V — R colors 
derived for the LSB component with the predic- 
tions of two different groups of evolutionary syn- 
thesis models: 

The Galaxy Isochrone Synthesis Spectral Evo- 
lutionary Library, GISSEL 96 (Leitherer et al. 
1996; Bruzual & Chariot 1993), which include sim- 
ple stellar populations (SSP) or IB with different 
metallicities from Z = 0.0004 to Z = 0.10. The 
age varies from to 20 Gyr, and Salpeter and Scalo 
IMFS are considered. We have also used the pre- 
vious version of the models (GISSEL 95) in order 
to include a constant star forming rate (CSFR) 
— only available for solar metallicity. 

The Spectrophotometric Population Synthesis 
Library for Old Stellar Systems (Vazdekis et al. 
1996), a set of stellar population synthesis mod- 
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els designed to study old stellar systems. They 
synthesize SSP or follow the galaxy through its 
evolution from an initial gas cloud to the present 
time, assuming an analytical functional form for 
the SFR (see Section 3.1 in Vazdekis et al. 1996), 
and including the chemical evolution. Two differ- 
ent IMFs are considered: a unimodal IMF, with a 
power law form whose slope is a free parameter, 
and a bimodal IMF, equal to the unimodal IMF 
for stars with masses above 0.6 M©, but with a 
reduced influence of the stars with lower masses. 

When comparing with the predictions of the IB 
approximation, both models give ages of ~ 5 Gyr 
for the older stellar component of Mrk 370. 

However, when changing the star forming his- 
tory, we obtained different results, the B — V and 
V — R colors could not be reproduced by GIS- 
SEL95 with a CSFR. Next, we applied Vazdekis 
et al. (1996) models, which include chemical evolu- 
tion. We considered Salpeter IMF, and the differ- 
ent v (y is a constant fixing the timescale of star 
formation, see Sect 3.1 in Vazdekis et al. 1996), 
provided by the models. We found the following 
results: when v = 1 the observed colors are incon- 
sistent with the predictions of the models; with 
v = 5 colors agree with ages between 9 and 11 
Gyr; and in the rest of the cases [y = 10, 20 and 
50) colors are consistent with ages between 4 and 
7 Gyrs. 

Optical colors are not very sensitive to the evo- 
lution of the older stars; for instance, B — V in- 
crease by w 0.1 mag when ages varies from 5 to 
17 Gyr and the same applies to V — R. Such vari- 
ations are of the order of the uncertainties in our 
measured colors. Although no firm constraints 
could be put on the age of the underlying stel- 
lar population, the above comparison with mod- 
els point to a minimum age of 5 Gyrs; but it is 
clear that in order to better constrain the age of 
the LSB host, it is indispensable to include near- 
infrared colors in the analysis. 

6. Discussion 

It turns out from the above results that the kind 
of spectrophotometric study we have carried out 
is fundamental in order to disentangle the stellar 
populations in BCDs and derive their star forming 
histories. However, as we have already pointed 
out in the introduction, most of the work carried 



out so far has focused on statistical analysis of 
BCDs samples, while only a few papers have been 
devoted to examining in detail the characteristic 
of individual objects: 

The two extremely compact and low-metallicity 
BCDs, SBS 0335-052 and Tololo 65, both classified 
as iO (LT86), have been the subject of recent spec- 
trophotometric studies by Papaderos et al. (1998, 
1999). They found in both cases that the prop- 
erties of the underlying component are consistent 
with those of a stellar population not older than 
one Gyr. 

Noeske et al. (2000) report on results from a 
spectrophotometric analysis of two iI,C (" Cometary 
BCDs" , LT86) Mrk 59 and Mrk 71. Spectral popu- 
lation synthesis models, in combination with color 
magnitudes diagrams and color profiles, yield a 
most probable formation age of « 2 Gyr for the 
old stars in both galaxies. Two other cometary 
BCDs have been studied by Fricke et al. (2001) 
- Tololo 1214-277 — and Guseva et al. (2001) - 
SBS0940+544. In the latter case, no compelling 
evidences favor either a young or an old age. 

Despite the fact that nE/iEs are the most com- 
mon type of BCD, there are very few studies, with 
the same quality, of objects belonging to this class. 
Steel et al. (1996) and Mendez et al. (1999) have 
performed spectrophotometric observations of the 
iE BCDs Haro 3 and II Zw 33 respectively. Al- 
though in both cases an underlying component of 
older stars was detected, their data did not per- 
mit to determine its properties, and both stud- 
ies limited the analysis to the young starburst. 
The ages and star-forming scenarios they found 
for these two objects are consistent with our find- 
ings. They also found metallicities comparable to 
the metallicity we have derived for Mrk 370. More 
recently, a comprehensive analysis of the iE BCD 
galaxy Mrk 86 has been published by Gil de Paz 
et al. (2000) and Gil de Paz, Zamorano, & Gallego 
(2000). They studied the properties of 46 individ- 
ual star-forming knots as well as the underlying 
population of older stars. By applying evolution- 
ary synthesis models, they found that three well 
defined stellar populations are present in Mrk 86: 
the star-forming regions, with ages between 5 and 
13 Myr, and no significant age or metallicity gra- 
dient, a central starburst with an intermediate age 
of 30 Myr, and an underlying population of stars, 
with a surface brightness profile that can be de- 
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scribed by an exponential, with no color gradients 
and an age between 5 and 13 Gyrs. 

In all the above studies the star-forming knots 
present a similar range of ages, consistent with our 
findings for Mrk 370. However, the galaxies ap- 
pear to be different with respect to the properties 
of their underlying population of stars: while iO 
and cometary BCDs seem to be less evolved sys- 
tems, the structural properties and ages derived 
for the LSB of Mrk 370 are quite similar to those 
presented by the iE BCD Mrk 86. This fact may 
imply that BCDs are objects which have in com- 
mon the presence of an active starburst, which 
dominates their optical properties but, when the 
characteristics of the underlying stellar component 
are taken into account, they form a more heteroge- 
neous class. This picture could naturally explain 
the different evolutionary stages derived for the 
galaxies. 

We argue that considering BCD as a class might 
be then inappropriate. Therefore, looking for ab- 
solute unifying scenarios — for instance, trying to 
link the whole BCD class with the different types 
of dwarf galaxies (when only a fraction of them 
could be actually connected) — may be specula- 
tive at this time. 

On the other hand, Noeske et al. (2000) sug- 
gested that an evolutionary sequence might ex- 
ist connecting the different BCD subtypes, the 
iI,C galaxies being a possible link between the ex- 
tremely young galaxy candidates (iO class) and 
the more evolved iE/nE BCDs. 

The information available at the present time 
does not permit to discard any of these hypothesis. 
Comprehensive studies of larger samples of BCDs 
are required in order to find out which scenario is 
the correct one. 

Our kinematic results suggest ordered motion 
around the north-south axis of Mrk 370, the kine- 
matic center being located close to the center of 
the external isophotes; this could imply that the 
gas is rotating coupled with the main body of 
the galaxy (the underlying population hosting the 
starburst). However, we should take into account 
that, because of our limited spectral resolution, 
more complex gas motions could escape detec- 
tion. In fact, Ostlin et al. (1999, 2001), analyzing 
high- resolution (velocity sampling ~ 5 km sec -1 ) 
Fabry-Perot observations of BCD galaxies, found 



that, in general, the velocity fields of the studied 
galaxies are irregular and distorted. Spectroscopic 
observations of Mrk 370 with high velocity resolu- 
tion are needed to better trace its gas kinematics. 

Besides, Mrk 370 present some other puzzling 
peculiarities: the peak in the continuum map does 
not coincide with the kinematic center, nor with 
the center of the outer isophotes, nor with one of 
the knots; and the position angle of the outer re- 
gions, seen in the i?-band grey-scale map, differs 
by about 30 degrees from the position angle dis- 
played by the Ha map. 

7. Summary and Conclusions 

Broad- and narrow- band images of the BCD 
Mrk 370, together with long-slit and two-dimensional 
spectra, have been analyzed in order to derive 
the properties of the different components of the 
galaxy and to constrain its evolutionary status. 
Our results can be summarized as follow: 

• Two different stellar components are clearly 
distinguished in the galaxy: the present star- 
burst and an underlying older population. 
The current star formation activity takes 
place in numerous knots, aligned in a north- 
east south-west direction. An extended LSB 
component, with regular appearance, hosts 
this star-forming area. 

• The oxygen abundance derived from the 
spectra is Z ~ Zq/2, a value relatively high 
for the BCD class; the high blue continuum, 
the absorption wings in the Balmer lines and 
the pronounced Balmer discontinuity in its 
spectrum are clear indicators of the presence 
of an evolved stellar population. 

• The starburst region is resolved into small 
and compact star-forming regions. We iden- 
tified a total of sixteen knots, and derived 
their individual properties. The brightest 
knots have luminosities similar to those of 
giant H II regions in M33 and in the LMC. 

The photometric parameters of the individ- 
ual star-forming knots were corrected from 
the contribution of emission lines, internal 
extinction and emission from the underlying 
LSB component. The colors before and af- 
ter the corrections are significantly different. 
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In Mrk 370 the internal extinction and the 
contribution of older stars are more impor- 
tant than emission lines in shifting the ob- 
served broad-band colors of the knots The 
color shifts are function of the position of 
the knot in the galaxy, and we strongly cau- 
tion against interpreting observed color gra- 
dients as age gradients before applying such 
corrections. 

We found that we can reproduce the colors 
of all the knots with an IB of star forma- 
tion and the Salpeter initial mass function 
(imf) with an upper mass limit of 100 M Q . 
The ages of the knots range between 3 and 
6 Myrs. 

• The age of the underlying stellar population 
has been estimated by comparing its optical 
colors with two different groups of evolution- 
ary synthesis models. Colors are consistent 
with ages larger than 5 Gyr. 

• The inner, brightest star forming-regions, 
previously cataloged as double nucleus, are 
normal H II regions. We would like to 
stress that the "multiple nuclei" often seen in 
broad-band images of BCD galaxies may be 
just the superposition of off-center luminous 
H II regions. If Mrk 370 is not a double- 
nucleus galaxy, then there may be no com- 
pelling evidence that it has experienced re- 
cent interactions. Indeed, there are no other 
signs of interactions; the regularity of the 
velocity field, and the fact that the central 
knots share the same kinematics of the main 
body point against a recent merger or in- 
teraction event. Also, Nocskc ct al. (2001) 
report that Mrk 370 is an isolated galaxy. 

• We obtained the mean ionized gas velocity 
field in the central part of the galaxy. The 
velocity field seems to be regular and sim- 
ilar to that of a rotating disk (west region 
approaching, east region receding) . 

All in all, spectrophotometric studies of BCDs 
are of paramount importance in order to derive 
their evolutionary status. These analyses must in- 
clude, in addition to the whole optical dataset, 
near-infrared photometry (NIR) — colors in the 
NIR better trace the properties of the older stellar 



component — as well as deep optical spectropho- 
tometry of the host galaxy. Besides, high reso- 
lution spectroscopy of BCDs will help trace the 
kinematics of the galaxy and provide valuable in- 
formation on the mechanisms triggering the star 
formation. 
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Table 1 
Basic data of Mrk 370 



Galaxy 


Other designations 


R.A. (1950) Dec. 


M B 


D (Mpc) 


Mhi 




Mrk 370 


NGC 1036, 
UGC 02160 


02 h 37 m 40 s 19° 05' 01" 


-17.20 


12.9 


0.36xl0 9 


4.2 x 10 9 



Note. — (1) Neutral hydrogen mass Mhi and total mass Mt in units of M Q ; both from Thuan and Martin 
(1981). 



Table 2 
Log of the observations 



Date 


Telescope 


Instrument 


Filter/grism 


Exposure time (s) 


Oct. 98 


JKT 1.0m 


Cass, focus 


U 


2400 


Oct. 98 


JKT 1.0m 


Cass, focus 


B 


1500 


Oct. 98 


JKT 1.0m 


Cass, focus 


V 


1000 


Oct. 98 


JKT 1.0m 


Cass, focus 


R 


800 


Oct. 98 


JKT 1.0m 


Cass, focus 


I 


1200 


Dec. 97 


CAHA 2.2m 


CAFOS 


6569 (113) 


5400 


Dec. 97 


CAHA 2.2m 


CAFOS 


6462 ( 98) 


5400 


Sep. 98 


INT 2.5m 


IDS 


R300V 


1800 


Aug. 99 


WHT 4.2m 


INTEGRAL 


600g/mm 


3600 



Note. — JKT = Jacobous Kapteyn Telescope, ORM (La Palma), Spain 
(detector: Tck 1024 x 1024, 0'.'33/pixel). CAHA = Ccntro Astronomico Hispano 
Aleman, Almeria, Spain (Site-ID 2048 x 2048, 0'.'53 /pixel) / INT = Isaac Newton 
Telescope, ORM (La Palma), Spain (EEV10 4100 x 2048 13. 5^ pixels). WHT 
= William Hcrschcl Telescope, ORM (La Palma), Spain (Tex 1024 x 1024 24[i 
pixels). 
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Table 3 

Reddening corrected line intensity ratios 



Line 


Ion 




A 






B 




Q 
O 




K A ) 






^X 


-W\ 


r\ 


-W x 




tx 


-Wx 


3727 


\0 III 


0.26 


2.74 ± 0.06 


35.1 ± 0.9 


3.22 ± 0.09 


136 ± 


10 


2.68 ±0.04 


48 ±1 


3869 


[Ne III] 


0.23 


0.50 ± 0.07 


5.1 ± 0.8 


1.1 ±0.2 


330 ± 166 


0.49 ±0.06 


7 ±1 


4101 


H<5 


0.18 


0.28 ± 0.03 


2.1 ± 0.4 


0.24 ±0.04 


7 ± 


1 


0.26 ±0.02 


2.5 ± 0.3 


4340 




0.14 


0.45 ± 0.02 


3.0 ± 0.2 


0.50 ±0.02 


14.0 ± 


0.8 


0.49 ±0.07 


5 ±1 


4861 


H/3 





1.00 ± 0.01 


12.6 ± 0.1 


1.00 ± 0.02 


26.9 ± 


0.8 


1.00 ±0.01 


16.0 ± 0.3 


4959 


[0 III] 


-0.02 


0.53 ±0.01 


6.1 ±0.2 


0.60 ±0.02 


15.5 ± 


0.6 


0.52 ±0.01 


7.6 ±0.1 


5007 


[0 III] 


-0.03 


1.60 ±0.02 


19.2 ±0.2 


1.93 ±0.03 


51.2 ± 


0.6 


1.63 ±0.02 


25.0 ±0.2 


5876 


Hel 


-0.23 


0.06 ±0.01 


0.8 ±0.1 


0.10 ±0.01 


3.5 ± 


0.3 


0.069 ±0.006 


1.3 ±0.1 


6548 


[N II] 


-0.34 


0.14 ±0.01 


2.5 ±0.2 


0.11 ±0.01 


4.2 ± 


0.6 


0.149 ±0.008 


3.4 ±0.2 


6563 


Ha 


—0.34 


2.85 ±0.05 


56.8 ±0.6 


2.93 ±0.09 


112 ± 


2 


2.86 ±0.03 


68.1 ±0.7 


6584 


[N II] 


-0.34 


0.50 ±0.02 


10.1 ±0.5 


0.49 ±0.02 


24 ± 


2 


0.476 ± 0.009 


11.6 ±0.3 


6678 


Hel 


-0.35 


0.04 ± 0.02 


0.8 ±0.4 


0.06 ±0.03 


3 ± 


1 


0.030 ± 0.008 


0.7 ±0.2 


6717 


[S II] 


-0.36 


0.50 ±0.01 


10.0 ±0.2 


0.45 ±0.02 


21.3 ± 


0.9 


0.471 ± 0.008 


11.5 ±0.2 


6731 


[S II] 


-0.36 


0.38 ±0.01 


7.4 ±0.3 


0.34 ±0.02 


17 ± 


2 


0.365 ±0.007 


9.0 ±0.2 


7135 


[Ar III] 


-0.41 


0.11 ±0.01 


2.4 ±0.4 


0.10 ±0.01 


4.9 ± 


0.8 


0.093 ± 0.006 


2.5 ±0.2 


C(H/3) 






0.31±0.02 


0.39±0.03 




0.285±0.004 


w abs (A) 






1.3 






0.0 




1.4 




F(H/3) 






2.67±0.0014 


1.55±0.12 




3.87 ±0.06 



Note. — Reddening-corrected line intensities (normalized to H/3 = 1) for the regions extracted from the long slit 
spectrum of Mrk 370. The integrated spectrum S is presented in the rightmost columns. Balmer lines are corrected 
from underlying stellar absorption. The reddening coefficient, C(H/3), the value of the absorption correction, W a b s , 
and the measured H/3 flux, F(H/3) (xl0~ 14 erg cm~ 2 sec -1 ) are also included. 



Table 4 

ha photometry of the individual knots in mrk 370 



Knot 


Flux(Ha) 


logi(Ha) 


log N Lym 


log(M H n/M Q ) 


-W(Ha) 


-iy(Ha) CO rr 


1 


5.70 


38.05 


49.92 


3.13 


201 


364 


2 


1.96 


37.59 


49.46 


2.67 


170 


534 


3 


5.08 


38.00 


49.87 


3.08 


138 


260 


4 


8.01 


38.20 


50.01 


3.28 


173 


329 


5 


4.13 


37.91 


49.78 


2.99 


45 


71 


6 


0.68 


37.13 


49.00 


2.20 


94 




7 


1.64 


37.51 


49.38 


2.59 


38 


116 


8 


37.30 


38.87 


50.74 


3.95 


58 


63 


9 


7.29 


38.16 


50.03 


3.24 


81 


119 


10 


4.56 


37.96 


49.82 


3.03 


114 


578 


11 


147.07 


39.47 


51.33 


4.64 


373 


533 


12 


232.17 


39.66 


51.53 


4.96 


192 


226 


13 


2.32 


37.66 


49.53 


2.74 


71 




14 


33.31 


38.82 


50.69 


3.90 


237 


579 


15 


4.01 


37.90 


49.77 


2.98 


73 


223 


16 


5.39 


38.03 


49.90 


3.11 


220 


656 



Note. — The data are corrected from interstellar extinction and from [N II] emission. The 
equivalent widths are shown before and after correction from the contribution of the underlying 
stellar emission to the underlying continuum (except knots 6 and 13). Ha fluxes arc in 10" 15 
erg cm" 2 sec" 1 units; Ha luminosities in erg sec" 1 units. VK(Ha) in Amstrongs 
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Table 5 

Broad-band photometry of the individual knots detected in Mrk 370. 



Knot 


D 
D 


U — n 


D — V 


V — H 


V — 1 


7~) 

n c 


IJJ D\ 

[U - B) c 


(B - V) c 


(V - R) c 


(V - J )c 


1 


20.19 


-0.57 


0.39 


0.25 


0.88 


19.67 


-0.94 


0.05 


0.04 


0.24 


2 


21.01 


-0.56 


0.33 


0.38 


1.01 


20.70 


-1.05 


-0.14 


0.23 


0.26 


3 


20.15 


-0.45 


0.47 


0.39 


1.00 


19.69 


-0.88 


0.12 


0.21 


0.47 


4 


19.88 


-0.54 


0.42 


0.44 


1.02 


19.43 


-0.97 


0.05 


0.30 


0.47 


5 


19.28 


-0.20 


0.53 


0.38 


0.96 


18.76 


-0.57 


0.24 


0.17 


0.45 


6 


21.82 


-0.23 


0.67 


0.31 


1.41 


22.24 


-1.54 


0.14 


-0.04 


1.20 


7 


20.11 


-0.09 


0.60 


0.41 


1.07 


19.88 


-0.74 


0.24 


0.17 


0.46 


8 


16.80 


-0.42 


0.28 


0.29 


0.66 


16.04 


-0.55 


0.06 


0.10 


0.33 


9 


19.17 


-0.20 


0.50 


0.40 


0.94 


18.56 


-0.47 


0.24 


0.18 


0.55 


10 


20.00 


-0.24 


0.57 


0.47 


1.13 


19.73 


-0.85 


0.21 


0.24 


0.66 


11 


17.46 


-0.65 


0.33 


0.39 


0.69 


17.37 


-0.79 


0.19 


0.28 


0.46 


12 


15.77 


-0.45 


0.31 


0.33 


0.66 


14.80 


-0.60 


0.02 


0.05 


0.21 


13 


20.13 


-0.11 


0.56 


0.48 


1.06 


19.97 


-0.84 


0.12 


0.10 


0.43 


14 


18.43 


-0.56 


0.35 


0.41 


0.82 


17.90 


-0.92 


-0.01 


0.14 


0.17 


15 


19.81 


-0.45 


0.42 


0.41 


1.25 


19.61 


-1.14 


-0.15 


0.26 


0.76 


16 


20.50 


-0.48 


0.32 


0.47 


0.84 


20.21 


-1.08 


-0.28 


0.45 


-1.68 



Note. — Columns 2 to 6: only the Galactic extinction correction has been applied; columns 7 to 11: values 
corrected from interstellar extinction, contribution of emission line and emission from the underlying host galaxy. 



Table 6 

Structural parameters characterizing the LSB host. 



Band 


Mo 


a 




M host 


-^SB/host 


B 


22.56 


1.06 


-16.69 


-16.13 


1.68 


V 


21.74 


1.05 


-16.91 


-16.93 


0.98 


R 


21.04 


0.96 


-17.32 


-17.44 


0.89 


I 


20.57 


1.14 


-17.23 


-18.28 


0.38 



Note. — Columns 2, 3: centrals surface brightness 
(mag arcsec~ 2 ) and scale-length (kpc); column 3, 4: 
absolute magnitude of the starburst and of the host 
LSB in each band; column 5: the ratio of the lumi- 
nosity of the starburst over the luminosity of the un- 
derlying component. As mentioned in the text, the 
{/-band image was too shallow to reach out the outer, 
starburst-free regions. 



Table 7 

Derived ages (in Myr) for the individual star-forming knots in Mrk 370 



Knot 


Agey-s 


A g e w(Ha) 


Knot 


Ageu-B 


A g e _Eiv(Ha) 


1 


3.2-4.8 


4.8-5.1 


9 


5.6-6.4 


5.9-6.3 


2 


3.1-3.8 


4.6-4.8 


10 


4.5-5.0 


4.6-4.8 


3 


3.5-5.0 


5.0-5.5 


11 


4.6-5.1 


4.6-4.8 


4 


3.2-4.7 


4.8-5.2 


12 


5.1-5.9 


5.1-5.7 


5 


5.1-6.0 


6.4-7.0 


13 


4.5-5.1 




6 




6.1-6.6 


14 


3.2-4.8 


4.6-4.8 


7 


4.8-5.3 


5.9-6.3 


15 


2.2-3.3 


5.1-5.7 


8 


5.3-6.2 


6.4-7.0 


16 


2.9-3.5 


4.3-4.7 
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Fig. 1— (left) Image of Mrk 370 through the Ha 
filter, taken with the 2.2m telescope at CAHA. 
The spatial distribution of the fibers (SB2) on the 
focal plane has been overlaid, (right) Fibers in 
the central array. Numbers indicate the actual 
position of the fibers at the slit. 
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Fig. 3. — Spectra of regions A and B, and the total 
integrated spectrum (region s). 



60 40 



40 -60 



a R.A. (arcsec) 



Fig. 2. — Contour plot of the continuum sub- 
tracted Ha image of Mrk 370. The slit position 
is indicated, and the two subregions selected in 
the long-slit spectrum are marked A and B. The 
individual star-forming regions are labeled. North 
is at the top and East is on the left. 
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Fig. 4. — [from left to right, top to bottom] (a) Mrk 370 continuum map obtained from the 2D spectroscopic 
data by integrating the signal in the indicated spectral interval; (b) V filter image; (c) Intensity map of the 
Ha emission line obtained by fitting a single gaussian to 2D spectroscopy data; (d) Ha filter image from the 
2.2m CAHA; (e) Two-dimensional distribution of [O III]/H/3 ratio; (f) Two-dimensional distribution of the 
extinction coefficient C(H/3). A cross marks the optical nucleus in each map. The geometrical center of the 
outer isophotes is indicated by a star. 
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800 




5 
Aq (arcsec) 



Fig. 5. — Velocity field of the ionized gas in the 
central region of Mrk 370, determined by averag- 
ing measurements from different emission lines. A 
black cross marks the optical nucleus. The kine- 
matic center is indicated by an X. The center of 
the outer isophotes is indicated by a black star. 




-*0 -ft 0. ?D. 




Fig. 7. — left panel: Continuum-subtracted Ha 
grey-scale map of Mrk 370; right panel: R-band 
image. North is at the top and East is on the left; 
axis units are arcsec. 
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Fig. 6. — Radial dependence of the systemic veloc- 
ity V S ys, amplitude of the velocity field projection 
on the galactic plane f2, and position angle of the 
major kinematical axis 5. Error bars only reflect 
the local rms in the measurements, and do not 
reflect the actual uncertainties, which are likely 
affected by systematic errors. 



a 




(U-B) 

Fig. 8. — Equivalent width of the Ha line ver- 
sus (U — B). We also plot the tracks of an in- 
stantaneous burst of z = 0.008, aged from 1 to 
30 Myr (Leitherer et al. 1999), for three IMFs: 
solid line, a=2.35, M up =100M o ; dotted line: 
a=3.30, M up =100M o ; short-dash line, a=2.35, 
M U p=30M Q . Triangles represent the observed 
(uncorrected) parameters of the individual knots; 
squares represent the knots parameters after ap- 
plying the corrections from reddening and from 
the contribution of emission lines and older stars. 
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0.4 0.2 0.2 0.4 0.6 OB 

(B-V) 

Fig. 9.— (U - B) versus (B-V). The track of 
an instantaneous burst of z = 0.008, aged from 1 
to 30 Myr (Lcitherer et al. 1999), are also plot- 
ted, for three different imfs: solid line, a— 2.35, 
M„ p =lOOM ; dotted line: a=3.30, M Mp =100M Q ; 
short-dash line, a=2.35, M up =3OM . Symbols 
are the same as in the previous figure 



8. APPENDIX A 



8.1. Atlas of Spectra 



The nuclear spectrum of Mrk 370 derived from 
the INTEGRAL data is shown in figure 10, in the 
whole wavelength range observed. We can easily 
recognize several emission lines. Here we only plot 
some of selected emission as spectrum diagrams 
(figures 11 through 14). Spectrum diagrams rep- 
resent line profiles from individual spectra at each 
point in the observed region and in a small spectral 
range. The spectra at each location are normal- 
ized by the local peak so as to better display the 
profile shape (lines nearer to the bright knots are 
brighter than those farther out). 




5500 6000 6500 

MA) 



Fig. 10. — Nuclear spectrum of Mrk 370, obtained 
by averaging the 7 fibers closest to the continuum 
maximum (r < 1".2). 
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Fig. 11. — Spectra distribution diagrams of the 
central 16" x 12" of Mrk 370. Two-dimensional 
distribution of the H/3 lines (the absorption wings 
are visible). The spectral range plotted is 4795- 
4950 A. 
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Fig. 12. — Two-dimensional distribution of 
[O III]A4959+[0 III]A5007 (4950-5045 A). 
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Fig. 14. — Two-dimensional distribution of 
[S II]AA6716,6731 (6720-6785) A). 
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Fig. 13. — Two-dimensional distribution of 
Ha+[N II]AA6548,6584 (6550-6640 A). 
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